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(57) Abstract 



This invention provides a signal processing method for en- 
hancing the communication quality and increasing the communi- 
cation capacity by reducing the effects of interference and noises 
with the nice beam pattern. The signal processing method pro- 
vides a beam pattern by computing an eigenvector corresponding 
to the maximum eigenvalue of an autocorrelation matrix of re- 
ceived signals in an antenna array system. The inventive signal 
processing method introduces a simplified computational tech- 
nique for generating a nice beam pattern having its maximum 
gain along the direction of the wanted signal and maintaining 
the gain toward the direction of the interfering signals in as low 
a level as possible. 
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A SIGNAL PROCESSING METHOD UTILIZING AN 
EIGENVECTOR CORRESPONDING TO THE MAXIMUM 
EIGENVALUE OF AN AUTOCORRELATION MATRIX OF 
RECEIVED SIGNALS FOR AN ANTENNA ARRAY SYSTEM 

Description 

Technical Field 

This invention relates to a signal processing technique for wireless 
communication systems, and more particularly to a signal processing 
method for reducing the effect of interference and noise by controlling 
beam patterns in real-time, in a communication system utilizing an 
antenna array. 

Background Art 

In general, an original signal transmitted by a certain transmitter 
(hereinafter, simply called a "wanted signal") is always received at a 
receiving set together with other plural interfering signals. Since the level 
of distortion in a telecommunication system is determined by the ratio 
between the power of the wanted signal and total power of all the 
interfering signals, even if the level of the wanted signal is much higher 
than each of the interfering signals, the distortion of the communication 
system can pose a serious problem when the total power of all the 
interfering signals proportionally increase according to when the number 
of the interfering signals is rather high. In conventional telecommunication 
systems, interfering signals make it very difficult to extract the 
information from the wanted signal. 

Although an antenna array system has been considered as a 
countermeasure to improve the problems caused by the interfering signals, 
no practical method of synthesizing the antenna array system in an actual 
telecommunication systems, particularly for mobile communication systems, 
has yet been suggested. The problems of applying conventional antenna 
array systems, which is based on the method of Eigen-Decomposition, is 
mainly due to its complexity and operating speed which is too large for 
real-time processing in telecommunication systems. 

The conventional technique about the antenna array system, was 
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introduced in the following references^ 

[1] M. Kaveh and A. J. Barabell, "The Statistical Performance of the 
MUSIC and Minimun-Norm Algorithms for Resolving Plane Waves in 
Noise," IEEE Trans., Acoust, speech and signal process., vol. ASSP-34, 
pp. 331-341, April 1986. 

[2] T. Denidni and G. Y. Delisle, "A Nonlinear Algorithm for Output 
Power Maximization of an Indoor Adaptive Phased Array," IEEE 
Electromagnetic Compatibility, vol 37, no. 2, pp. 201-209, May, 1995. 

[3] B. G. Agee, S. V. Schell, and W. A. Gardner, "Spectral 
Self-Coherence Restoral: A New Approach to Blind Adaptive Signal 
Extraction Using Antenna Arrays", Proc. of IEEE, Vol. 78, No. 4, pp. 
753-767, Apr 1990. 

The problems in most conventional methods of designing antenna 
array systems are, first, it, except the method introduced in [3], requires 
some knowledge about the location of the wanted signal apriori, and 
second, it requires so many computations that real-time processing cannot 
be performed. Especially, when the arrival angle of the wanted signal or 
the total number of signal sources is unknown, the required amount of 
computation becomes even larger, which makes it impossible to apply the 
conventional method of synthesizing the antenna array system to practical 
signal environments,, such as mobile communications. 

Another undesirable feature of most conventional methods of 
designing antenna array systems is that the performance and/or the 
complexity of the system to be built is affected by the coherence and/or 
cross correlation of the wanted signal with respect to the interfering 
signals. This means that the antenna array system often requires 
additional complexities when: the signals are fully or partially coherent so 
that the resultant system becomes too complicated for real-time 
processing, which is very important, especially in mobile communications. 
Another conventional methods, such as the one shown in [3], do not work 
at all if the wanted signal cannot be extracted from the interfering signals 
with a particular frequency separation at a proper time lag. 

This invention introduces a new signal processing technology of 
designing an antenna array system that provides for a nice beam pattern 
having its maximum gain along the direction of the wanted signal 
maintaining the gain along the directions of interfering signals in 
relatively much lower levels. Under the assumption that the wanted signal 
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is sufficiently larger in magnitude than each of the interfering signals, the 
proposed technique generates the desired beam pattern without requiring 
any knowledge about the wanted signal as well as the interfering ones. 
Another important and attractive aspect of the proposed technique is that 
the total required amount of computation is so small that the optimal 
parameters of the antenna array system are produced on a real-time 
basis. In fact, the signal processing apparatus, which forms the 
beamforming module of the antenna array system introduced . in this 
invention, can easily be implemented with a normal, off the shelf, digital 
signal processor. 

The primary objective of this invention is to introduce a new signal 
processing method for designing a beamforming module of ah antenna 
array system in order to apply it at the base station of a mobile 
communication system for receiving and transmitting the signal of each 
subscriber in a cell with a nice beam pattern which is provided 
individually for each subscriber of the cell. The proposed technique can 
also be applied in other signal environments such as WLL( wireless local 
loop) and other fixed communications as well as mobile communications. 

Disclosure of Invention 

To solve the above mentioned problems, it is an object of the present 
invention to provide a signal processing method for enhancing the 
communication quality and increasing the communication capacity by 
reducing the effects of interference and noises with the nice beam pattern. 
The signal processing method presented in this invention provides a beam 
pattern by computing an eignevector corresponding to the maximum 
eigenvalue of an autocorrelation matrix of received signals in an antenna 
array system. The signal processing method introduces a simplified 
computational technique for generating a nice beam pattern having its 
maximum gain along the direction of the. wanted signal and maintaining 
the gain toward the direction of the interfering signals in as low a level 
as possible. 

To accomplish the object of the present invention, there is disclosed 
a signal processing method, comprising the steps of: (a) setting the initial 
estimation vector m(0) with a predetermined value, setting the initial gain 
vector m>(0) with the normalized estimation vector, i.e., m(Q) = z(Q) and 
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w(0). = \\m(0)\\ ' and setting the snapshot index n with 1 T i.e., « = 1; (b) 
updating the autocorrelation matrix JZ x (n) with a new signal vector x(n) 
received at the present snapshoht by JS x (n) = fB x (n—\) +. (1—/) 
^U) z"(n), where the initial matrix is set with the initial signal vector 
by JS X (Q) - jc(0) x H (Q) and the forgetting factor / is properly set with a 
value between 0 and 1, i.e., 0 <: / <1; (c) updating the estimation vector 

M(n) and gain vector jsiin) from the new autocorrelation matrix by 

ji'(n) - JSx(n) Mt(n-l) and Mn) = \\%(*)\\ > respectively; (d) going to 

the following step (e) t for more iteration in the present snapshot, or going 
to said step (f), for no more iteration in the present snapshot; (e) going 
back to said step (c), with ju(m— ■ Mii'n) for more iteration in the 
present snapshot; and (f) generating the final array output y(n) at the 
present snapshot by y(n) = w H (n) x(n) , and going back to said step, 
(b), with incremented snapshot index «, i.e., n «— for continuing the 

procedure at the next snapshot. 

Also, to accomplish the object of the present invention, there is 
disclosed a signal processing method, comprising the steps of: 

(a) setting the initial estimation vector j>(0) with a predetermined 
value, setting the initial gain vector j^.(0) with the normalized estimation 

vector, i.e., ^(0) = jf(0> and mi(0) = 9 and setting the sna P shot 

index n with 1, i.e., w=l; 

(b) Receiving a new signal vector x{n) at the present snapshoht; 

(c) " updating the estimation vector j/(w) and gain vector win) from 
the new signal vector by jtin) = f jj(n-l) + (1—/) sin) z H (n) win — I) 

and M(n) = 1 res P ect i ve ^y» where the forgetting factor / is 

properly set with a value between 0 and 1, i.e., 0 <> / <1; 

.■■(d) going to the following step (e), for more iteration in the present 
snapshot, or going to said step (f), for no. more iteration in the present 
snapshot; 

(e) going back to said ste;p (c), with ji(n-l) and jyin-l) 
<- win) for more iteration in the present snapshot; 

(f) generating the final array output y( n) at the present snapshot by 
yin) = M H (n) \x(n) , and going back to said step (b), with incremented 
snapshot index n, i.e., n +- for continuing the procedure at the next 
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snapshot. 

Also, to accomplish the object of the present invention, there is 

disclosed a signal processing method, comprising the steps of: (a) setting 

th-2 initial estimation vector tf(0) with the initial signal vector ^(0) , 

setting the initial gain vector w(0) with the normalized estimation vector, 

i.e., jKO)- = ^(0) and j#(Q) = ||^(0) || ' anc * sett * n S the snapshot index « 

with 1, i.e., « = 1; (b) receiving a new signal vector &(n) at the present 
snapshoht; (c) updating the estimation vector v(n) and gain vector w(n) 
from the new signal vector by m(h) = / v(n-l) + (1— f) An) x H (n) 

w(n-\) and win) = "jf7(^|f . respectively, where the forgetting factor / 

is properly set with a value between 0 and 1, i.e., 0^/<l; (d) r generating 
the final array output y(n) at the present snapshot by y(n) - ju H (n) 
x(n) , and going back to the second step, (b), with incremented snapshot 
index n, i.e., n +- for continuing the procedure at the next snapshot. 

Brief Description of Drawings 

The novel features believed characteristic of the invention, as well as 
other features and advantages thereof, will best be understood by 
reference to the following detailed description of a particular embodiment, 
read in connection with the accompanying drawings, wherein: 

Fig. 1 is a block diagram of an antenna array system that adopts a 
signal processing method according to the present invention; - 

Fig. 2A is a flowchart that descibes the four steps of computing the 
weight vector m according to the first embodiment of the inventive 
method;. 

Fig. 2B is a flowchart that descibes the six steps of computing the 
weight vector m according to the second embodiment of the inventive 
method; 

Fig. 2C is a flowchart that descibes the four steps of computing the 
weight vector w according to the third embodiment of the inventive 
method; 

Fig. 2D is a flowchart that descibes the six steps of computing the 
weight vector jui according to the fourth embodiment of the inventive 
method; 
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Fig . 3A shows the number of clocks required to obtain the target 
weight vector by utilizing a microprocessor; and 

Fig. 3B shows the SIR and SNR performances of the antenna array 
system designed by the method presented in this invention. 

Best Mode for Carrying out the Invention 

A preferred embodiment of the present invention will be explained 
below with reference to the accompanying drawings. 

The signal processing method that is proposed in this invention 
generates a proper gain vector (which is often called weight vector as 
well), that provides a beam pattern having its maximum gain along the 
direction of the wanted signal maintaining the gain in the other directions 
in as low a level as possible. This can be accomplished in one of two 
approaches. 

The first approach for providing the nice beam pattern is to optimize 
the value of the complex gain that is to be multiplied to each signal 
received at each antenna element, and the other approach is to optimize 
the value of the phase delay that is to be added to each signal received 
at each antenna element. Since each element of the gaion vector in the 
first approach is to be weighted (multiplied) to each element of the signal 
vector, the gain vector is often referred to as "weight vector" as well. 

Only the first approach is explained in detail in this document because 
of the following two reasons^ first, those two approaches are 
mathematically equivalent, second, the second approach is more costy in 
hardware and the beamforming performance of the second approach is 
inferior to the first one in terms of angle accuracy. However, the second 
approach can easily be implemented by utilizing the procedure explained in 
this invention. Therefore a signal processing method that is implemented 
in the second approach should be a part of this invention as long as it 
utilizes the idea introduced in this document 

In other words, this invention determines the complex gain vector "w" 
in such a way that the desired beam pattern is formed, and as a result, 
the output of the array antenna system, i.e., the Euclidean inner product 
of the signals induced at the antenna elements and the complex gain 
vector, should be as close to the wanted value as possible. 
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If the magnitude of every element of the complex gain vector is 
normalized to l t to multiplying the signal received at each antenna 
element by the corresponding element of the complex gain vector w is 
equivalent to adding the phase delay to the signal by the amount of the 
phase term of each corresponding element of the complex gain vector. 
Therefore, multiplying the signal vector by the gain vector is equivalent < 
to adding the phase of the signal vector by the amount of the phase term 
of the gain vector. 

As mentioned above, the same effect can also be obtained by 
appending the time delay to the signal received at the i_th antenna 
element by the amount of 0i divided by 2tt f c , where <f>\ and f c denote 
th? phase delay to be added to the signal received at the i_th antenna 
element and the carrier frequency, respectively. 

For a linear array system having a uniform spacing of r4jr- between 

adjacent antenna elements, where Ac denotes the wavelength at the 
carrier frequency, the signal induced at the m_th antenna element can be 
represented after the frequency down conversion as follows:. 

x m (t) = k |: i S k (t)e ,(m - 1) ' sin ^ + n m (t) (1) 

where 6 k denotes the incident angle of the k_th signal and Sk(t) is 
the k_th transmitted signal observed at the receiving end. 

The subscript m in equation (1) represents the antenna element. The 
reference antenna element is assigned to be m=l and the other antenna 
elements are assigned the next numbers, i.e., m=2, 3, in the order of 
the magnitude of the phase of the signal induced at each antenna element. 

In eq. (1), one of the M - signals is the wanted signal. For example, 
when Si(t) is the wanted signal it must be received at the antenna array 
system while all the other M-l signals, i.e., S2(t), S3(t), Sm(0, are 
interfering signals to be rejected, together with the noise n m (t), for a good 
signal reception. Although eq. (1) is valid for the linear array with the 
uniform half-wavelength spacing, the technique provided in this invention 
can be generally applied to non-uniform spacing or non-linear array 
systems as well. 

For non-uniform spacing arrays, if the distance of the m_th antenna 
element from the reference antenna element is d m , then there exists a 



7 



WO 98/00930 



PCT/KR97/00129 



phase difference in the signal induced at the m_th antenna element by 
2 ^-^- sin <9 k compared to the phase of the signal at the reference 

antenna element. Thus, the signal induced at the m_th antenna element 
for non- uniform and/or non-linear array systems can be written as 
follows: - 

x m (t) - k 2 i S k (t) e c + n m (t). 

In this invention, in order to make the phase delay to be appended at 
each antenna element be a positive quantity, the reference antenna element 
is defined as the antenna element at which the induced signal has the 
latest phase in the receiving array. In the transmitting array system, 
therefore, the antenna element at which the induced signal has the earliest 
phase is the reference antenna element. Defining the reference antenna 
element in the way explained above, the array antenna system can easily 
be designed by appending the zero phase delay to the signal at the 
reference antenna element and the proper positive amount of the phase 
delay to the signal at the other antenna elements. For an array antenna 
system consisting of N antenna elements, the array receives the N-by-1 
signal vector at every snapshot The autocorrelation matrix of the received 
signals can be written as shown in eq. (2). 

The term "snapshot" in this document denotes the time period during 
which the new gain vector (or, phase delay vector) is computed upon 
receiving the new signal vector. In this invention, the array antenna 
system that adapts to the new signal vector can be designed at each 
snapshot by determining the proper gain vector (or, phase delay vector) 
for each new signal vector received at every snapshot. 

Rx(J) = j ^-^a + lTs) _x H (t + lT s ) (2) 

where the underlined quantities denote the vector or matrix, Ts is the 
snapshot period and superscript H is the Hermitian operator. The N-by-1 
signal vector x(t), of which the number of elements is N, consists of the 
received signal x m (t) for m. = 1, 2, N, which is explained in eq.(l) as 

follows : 

. ^(0 = [x\(t)x 2 (t) - x N (t)] T (3) 
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where superscript T denotes the transpose operator. 

However, eq. (2) is valid only when the arrival angles of all the 
signal components remain unchanged. In a time-varying environment 
where each signal source moves during the communication, as in the 
mobile communication environment, the autocorrelation matrix cannot be 
obtained by eq. (2) because the arrival angle of the signal source changes 
at every snapshot. 

Therefore, in time-varying signal environments, it is recommended 
that the autocorrelation matrix be computed in an iterative manner as 
follows* 

R X (J + 1) = f -Rx(J) + ^((J + 1)T S ) _x H ((J + l)t s ) (4) 
where R X (J + D and E X (J) denote the autocorrelation matrix at the 

J+lst and J_th snapshot, respectively, and f denotes the forgetting factor 
in the range between 0 and 1. Since communication environments, 
especially mobile communications, are generally time-varying 
environments, the autocorrelation matrix in this invention is computed by 
eq. (4) rather than eq. (2). 

From various computer simulations, it is recommended to set the 
value for the forgetting factor, /, in the range between 0.8 and 0.99 for 
optimal performances in land mobile communications. 

Now, the design of the optimal array antenna system will be 
explained in more detail by taking the practical examples of actual 
applications. The eigenvalues of the autocorrelation matrix, determined 

by eq. (2) or (4), can be sorted by the magnitude as A i^ A 2^ - - • > A 
n. The largest eigenvalue A 1 is determined by the signal components, not 
the noise components, regardless of the number of signal sources or 
antenna elements. 

Therefore, the eigenvector corresponding to the largest eigenvalue A \ 
exists in the signal subspace as follows: 

_£ = f\£ i^(^i) (5). . 

where the complex quantity ' $ i is a constant determined by the 

magnitudes and distribution of the wanted and interfering signals, and the 
vector „a( e i) is the steering vector of the i_th signal component in the 
following form: 

9 
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/ n -\ r i j n sin 6 i j(N- 1) k sin e -I /C v 

_a( 0 1") = LI e • • • e J. (6) 

Now, suppose the magnitude of the wanted signal is sufficiently 
larger than each of the interfering signals such that the condition shown 
in (7) is satisfied. ^ 

I Si(t) | > > I Si(t) | for i*=l (7) 

In a signal environment in which condition (7) is satisfied, the 

eigenvector A \ corresponding to the largest eigenvalue ,can be 
approximated as: 

jsli = <? i Jk(0 i)> (8) 

This means that the steering vector, a( 6 x ), of the wanted signal is 

almost the same as the eigenvector corresponding to the largest 
eigenvalue, except that the complex- valued constant, $ i, is multiplied. 
Therefore, under the condition that the wanted signal is sufficiently larger 
than each of interfering signals, the maximum gain of the array antenna 
system will approximately point to the direction of the source of the 
wanted signal if the gain vector to be appended to the antenna elements 
of the array system is determined by the eigenvector corresponding to the 
largest eigenvalue of the autocorrelation matrix of the signals impinging 
upon the array system. 

In conclusion of the above discussions, this invention suggests that 
the gain vector can be determined by the following equation: 

Now, the practical way of computing the optimal weight vector is 
presented. As mentioned previously, under a particular signal environment 
where the wanted signal is sufficiently larger than each of interfering 
signals, the array antenna system having the desired beam pattern, which 
provides the maximum gain along the direction of the wanted signal 
source, can be obtained by taking the weight vector w with a normalized 
eigenvector _£! corresponding to the largest eigenvalue A \ of the 

autocorrelation matrix. However, to obtain the autocorrelation matrix itself 
requires a lot of computations, as shown in eqs.(2) and (4). Moreover, it 
is not a simple task to compute the eigenvector corresponding to the 
largest eigenvalue of the matrix. What makes the problem even more 
complicated is that the arrival angle of each signal changes at every 
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snapshot in mobile communications such that the eigenvector to be 
obtained varies at every snapshot. 

Considering the time-varying signal environment, this invention 
introduces a method of computing the weight vector w with the 
approximated value for the eigenvector jz x in. an iterative way. This 
means that the weight vector obtained at the last snapshot is updated at 
the next snapshot upon the reception of a new signal vector that contains 
time- varying statistics mainly due to the different arrival angles at each 
snapshot. The iterative procedure for determining the weight vector with a 
normalized eigenvector corresponding to the maximum eigenvalue of the 
autocorrelation matrix is described in detail in the remainder of this 
section. 

In general, the eigenvector corresponding to the maximum eigenvalue 
of a matrix can be obtained by the method known as "Power Series 
Method" (abbreviated as "PSM" in the remainder of this document), of 
which the basic concept can be found in many textbooks of linear 
algebra. Since the signal environment considered in this invention is such 
that the wanted signal is much larger than each of interfering signals, the 
above mentioned technique of PSM can produce the target eigenvector 
with relative efficiency. The mathematical details about the PSM is given 
in the following textbook: 

[4] Howard Anton, Elementary Linear Algebra- Chapter 8, Section 3, 
John Wiley & Sons Inc., 1984. 

The key part of this invention is to provide a simple procedure for 
computing an approximated eigenvector of the largest eigenvalue of the 
autocorrelation matrix that must be updated at each snapshot upon the 
reception of a new signal vector of the antenna array system. For a given 
matrix £ t the procedure of determining the weight vector with the target 
eigenvector by means of the PSM can be summarized as given in [4] as 
follows: 

<step 1> Set an initial estimation vector i<0) avoiding zero vector. 
<step 2> Update the estimation vector by multiplying the present 
estimation vector An) to the autocorrelation matrix, i.e., 
v(n+ 1) = £ An) . 

<step 3> Determine whether or not the computation of the estimation 
vector is repeated for more accuracy. Go back to <step 2> with An) «— 
An+1) for repeatition, otherwise, go to the following step, i.e., <step 4>. 

n 
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<step 4> Generate the weight vector j£>(«+l) by normalizing the 
estimation vector obtained in the last step, i.e., 1) - + 

\h>(n+l)\\ . 

The eigenvector corresponding to the largest eigenvalue can be 
obtained in a finite number of iterations from the above- written procedure 
unless there exists more than one larger eigenvalues in the given matrix. 
Also, as mentioned previously, since the wanted signal is assumed to be 
much larger in magnitude than each of interfering signals, the target 
vector can be obtained in relatively small number of repeatitions., A 
serious problem, however, in the above- written procedure is that a matrix 
multiplication is required in computing the target eigenvector, as shown in 
<step 2>. It means that, if L_times repeated computations for a given 
matrix is required, then the amount of computation due to the matrix 
multiplication in <step 2> would be proportional to 0(LN ). Furthermore, 
the above- written procedure consisting of 4 steps is applicable for 
obtaining the target eigenvector only after the autocorrelation matrix has 
been computed. Therefore, in order to apply the PSM to compute the 
weight vector, some additional steps must be performed successfully for 
receiving a signal vector and updating the autocorrelation matrix based on 
the new signal vector prior to the above-written procedure of PSM. 
Consequently, it seems to be impossible to apply the procedure of PSM, 
shown above, to the problem of designing a practical antenna array 
system, considering the total required amount of computation for the 
reception of a new signal vector, computation of an autocorrelation matrix 
based on this new signal vector, and performing the above-mentioned 4 
steps of PSM. In order to compensate for the disadvantages of the 
original PSM mentioned above, this invention provides a simplified 
procedure in which the matrix operation in <step 2> can be approximated 
with proper vector operations. Also, the method proposed in this invention 
suggests an efficient way of setting the initial value of the estimation 
vector ^>(0) from the signal vector received at the initial snapshot ^(0). 

As can be verified easily, the required number of iterations for 
computing the target eigenvector for a given matrix can significantly be 
reduced by selecting the signal vector itself as the initial value for the 
estimation vector at the very first snapshot, i.e., _fc<0) = x(0) . The basic 
reason is that the signal vector ^(0) itself must be a good approximation 
for the eigenvector corresponding to the largest eigenvalue at the first 
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snapshot because there must only be one signal eigenvector, meaning that 
the rank of the autocorrelation matrix at the very first snapshot is one. 
Starting from the initial estimation vector being set with the signal vector 
itself, i.e., AO) = .^(0) , the required number of iterations at each snapshot 
for achieving a given accuracy can be significantly reduced. In fact, it has 
been found in various computer simulations that the target eigenvector 
can be obtained with a reasonably good accuracy in only one "iteration at 
each snapshot, which means no repeated computation is required at each 
snapshot once the initial value for the stimation vector is set with the 
signal vector at the initial snapshot. 

Another defective feature of the original procedure of the PSM given 
in T4] is that, as mentioned previously, the matrix multiplication must be 
performed at step 2 at each iteration for a given matrix. Since the matrix 
itself must be updated upon the reception of a new signal vector, the total 
required amount of computation seems to be too large for real-time 
processing of practical mobile communication signals. 

The signal processing method for computing the eigenvector 
corresponding to the largest eigenvalue proposed in this invention 
improves the above-mentioned procedure of the PSM in such a way that 
the target vector is obtained in an adaptive manner with a relatively 
much simpler computational load. 

The autocorrelation matrix is updated whenever a new signal vector 
is received, as shown in eq.(4) which is rewritten as- 

£ x (n) = f £An-l) + (1-y) *{n) V(n) (4) 

where /.is a forgetting factor of which the value is predetermined in 

between 0 and 1, i.e., 0 / < L The new estimation vector An) is 

obtained by multiplying the weight vector computed in the last previous 

snapshot _u>(n-l) by the matrix updated with a new signal vector as 

shown in eq (4), i.e., j&n) = & x (n) Mn-\) . If the weight vector of the 

previous snapshot was computed successfully, since the magnitude of the 

wanted signal is much larger than each of interfering signals in the signal 

environment considered in this invention, the weight vector Mi(n) 

approaches the eigenvector corresponding to the largest eigenvalue of the 

autocorrelation matrix through the repeated procedure of the matrix 

multiplication ji(n) = _B x (n) jAn-Y) and normalization of-, the estimation 

vector win) = jAn)/\jAn)\. The required number of repeatitions of matrix 

multiplication and normalization depends on the accuracy of the initial 
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setting for the estimation vector jy(0) and how much larger the maximum 
eigenvalue compares to each of the other eigenvalues. As mentioned 
earlier, this invention suggests that the initial value for the estimation 
vector be set with the initial signal vector, i.e., m(0) = jc(0) . 

The above statements present a procedure of computing the weight 
vector to be applied as a gain vector to antenna elements in an adaptive 
manner. This means that the weight vector is obtained at each snapshot 
based on the solution computed in the last snapshot starting from the 
initial setting at the very first snapshot. However, as" for the complexity 
of the entire procedure, it still includes the matrix operation jy(«> = £ x (n) 
jt'vw-l), which should be performed as many times as the number of 
iterations at each snapshot. 

Developing a simplified procedure considering mainly the matrix 
operation An) "= J£ x in) win- 1) , the matrix operation can be written as 
foMows: 

An) = £ x {n) Mn-\)= (/ & x (n-l) + (l-/> i(n). & H { n) ) jd n ~ 1) . (10) 

in the meantime, in computing the autocorrelation matrix at each 
snapshot, as in eq.(4), it has been found in various computer simulations 
that the term (1-f) can be ignored without affecting the performance of 
the system. Thus, in the remainder of this document, the term (1-/) will 
net be shown in equations. 

The two terms of the right side of eq.(10), i.e., / JS,(«-1) w(n- 1) 
and z(n) £ H (n) w{n~ 1) , can be thought separately as follows. 

If the weight vector at the last previous snapshot _&;(« — 1) was 
computed successfully, the first term R x (n-\) jdn-1) can be written as 
JBx(«-D Mn-l) ^ A, Mn-l) (11) 

where is the maximum eigenvalue of the autocorrelation matrix at 
the previous snapshot £ x (n-l) \ But, once the weight vector win- 1) 
approximates to an eigenvector corresponding to the maximum eigenvalue 
of the matrix at the previous snapshot, A x is approximately equal to 
l|j<«- 1)11. Then, eq.(ll) can be written as 

£An-\) Mn-l) A\ Mn-l) * -|Lb(w— 1)11 jcK» — 1) = Jl(n-1). (12) 

Therefore, eq.(lO) that includes a matrix multiplication can be 
simplified into a vector equation as follows: 

Ji(n) = EAn) ui(n-\) * f A n- 1) + x(n) z H {n) _a<»-l). (13) 

The procedure of computing a normalized eigenvector corresponding to 
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the largest eigenvalue of the autocorrelation matrix at each snapshot can 
be summarized in a simplified form as follows: 

<step 1> Set the initial estimation vector jKO) and initial gain vector 
ji»(0) with the initial signal vector ^(0) and the normalized initial 

estimation vector, respectively, i.e., j;(0) = ^(0) and _a>(0) = \\^(q)\\ ' 

and set the snapshot index n with 1, i.e., n-\. 

<step 2> Receive a new signal vector at the present snapshot. 

<step 3> Update the estimation vector jiin) and gain vector Mi(n) 
from the new signal vector by jt(n)' =. J ji(n-l) + (\-f) x(n) x u (n) 
w(n-l) and ju(n) = 9 res P ectivel y' where the forgetting factor / 

is properly set with a value between 0 and 1, i.e., 0 ^ /<1. 

<step 4> Go to <step 5> for more iteration in the present snapshot, 
or go to <step 6> for no more iteration in the present snapshot. 

<step 5> Go back to < step 3> with jKw-1) *- ji(n) and j^(n-l) 
*- win) for more iteration in the present snapshot. 

<step 6> Generate the final array output y(n) at the present snapshot 
by y(n) = x(n) , and go back to <step 2> with incremented 

snapshot index n, i.e., »<- for continuing the procedure at the next 

snapshot. 

The repeated iteration at a snapshot may not be required in some 
signal environments, e.g., when the wanted signal is sufficiently larger 
than each of interfering signals. In this signal environment, the above 
procedure can be set more simply as follows^ 

<step 1> set the initial estimation vector _i/(0) and initial gain vector 
jjf(0) with the initial signal vector ^(6) and the normalized estimation 

vector, i.e., = jc(0) and _aK0) = \\%[q)\\ > and set the snapshot 

index n with 1, i.e., n-1. 

<step 2> Receive a new signal vector j;(n) at the present snapshot. 

<step 3> Update the estimation vector and gain vector j^Cn) 

from the new signal vector by ji(n) =. / ji(n-l) + (l-f) x(n) z H (n) 

win- 1) and = > respectively, where the forgetting factor / 

is properly set with a value between 0 and 1, i.e., 0^/<L 

<step 4> Generate the final array output y(n) at the present snapshot 
by y(n) = w H {n) jc(n) , and go back to <step 2> with incremented 
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snapshot index n, i.e., » «- «+l f for continuing the procedure at the next 
snapshot. 

It has been found in various computer simulations that the 
signal-to-(total)interference ratio (SIR) of the antenna array system 
adopting the proposed PSM can be improved by up to 8 or 10 dB 
compared to a single antenna system when the number of antenna 
elements in the array system is 12. The major factors that determine the 
amount of improvement in SIR are, firstly, whether the weight vector is 
produced from the autocorrelation matrix obtained by eq.(4) at every 
snapshot or it is produced through the simplified procedure shown in 
eq,(13), and secondly, whether or not there are enough number of 
iterations for computing the weight vector is assigned at each snapshot. 
The signal -to^ noise ratio (SNR) is improved as the number of antenna 
elements of the array system increases, i.e., the SNR at the output of the 
array system is N times the input SNR. 

Fig. 3 illustrates the performance of the antenna array system 
designed by the PSM presented in this invention. Fig. 3A shows the 
number of clocks required to obtain the target weight vector by utilizing a 
microprocessor. In the simulations for the data of Fig. 3A, a digital signal 
processor,for example TMS320C30, has been selected arbitrarily. The 
required number of clocks in Fig 3A is shown in comparison to another 
technique, namely a modified conjugate gradient method, that has been 
developed most recently. Fig. 3B shows the SIR and SNR performances of 
the antenna array system designed by the method presented in this 
invention. As illustrated in the figure, the proposed system provides about 
8dB improvement in SIR. The SNR, as mentioned earlier, is improved 
linearly as the number of antenna elements in the array increases. From 
the performance shown in Fig 3, it! is apparent that the antenna array 
system designed by the proposed method is excellent in terms of both 
performance and complexity. 

The amount of computation required to obtain the weight vector by 
the proposed method is, as can be counted from the procedure shown 
previously, about CK2N 2 +2N) . or Q(5N), depending on whether the 
autocorrelation matrix is actually computed by eq (4) or the approximation 
procedure of eq.(13) is invoked at each snapshot, respectively. As 
mentioned previously, the repeated computation of the weight vector at 
each snapshot is not necessary when the initial estimation vector is set 
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with the initial signal vector, as suggested in this invention. The primary 
reason for the repeated iteration at a snapshot not to be necessary is that, 
as mentioned previously, the wanted signal is assumed to be much larger 
than each of interfering signals. This condition is naturally satisfied in 
CDMA (Code Division Multiple Access) mobile communications because, 
at each receiving card of the base station, the received signal/consisting 
of every subscriber's transmitted signal, is cross-correlated with the! code 
sequence of which the subscriber's, signal is to be received. Therefore, 
af:er performing the cross -correlation with the code sequence of a 
to-be-received subscriber, the wanted signal must be much larger than 
each of the other subscribers' signals. If the proposed method of 
designing the antenna array system is applied to CDMA base station, the 
required condition is automatically satisfied at the output of the 
cross -correlator of each receiving card. More exactly speaking, in CDMA 
mobile communications, after the cross-correlation, the wanted signal 
becomes larger than each of interfering signals by PG times, where PG 
denotes the processing gain of the CDMA system, i.e., the ratio between 
the chip rate of the code sequence and bit rate of the original data. 

In order to implement a total antenna array system that encounters 
bcth the receiving and transmitting modes, the weight vector computed 
during the receiving mode by utilizing the procedure shown in this 
invention can be used during the transmitting mode. As mentioned earlier, 
the signal processing method of computing the weight vector of an 
antenna array system presented in this invention can significantly enhance 
the communication quality and increase the communication capacity. It can 
also extend the life time of the battery at each subscriber, because the 
base station adopting the proposed antenna array system receives the 
signal from each subscriber very efficiently with the nice beam pattern 
which provides its maximum gain along the direction of each target 
subscriber. 

Now, a more detailed explanation of the proposed apparatus and 
method will follow by taking practical examples. 

Fig. 1 illustrates a conceptual block diagram of an antenna array 
system that adopts a signal processing method disclosed in this invention. 
As shown in the figure, the antenna array system according to the 
present invention comprises an antenna array 1 consisting of plural 
antenna elements 11, a signal-receiving apparatus 7 generating said signal 
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vector z at each snapshot, and a signal processing apparatus 8 generating 
said weight (or, gain) vector w. ' 

The antenna array 1 consists of plural^ antenna elements. The signal 
induced at each antenna element encountering a phase delay in accordance 
with the incident angles includes all the signals each of which is 
transmitted from a subscriber in the cell. The signal-receiving apparatus 
7 takes the signals from the antenna array 1. The signal-receiving 
apparatus 7 includes common signal-receiving procedures such as low 
ncise amplification, frequency-down conversion, quasi-quadrature detection, 
etc, that are necessary to generate the signal vector x, each element of 
which represents the in-phase and quadrature component of the received 
signal. When the method proposed in this invention is applied in CDMA 
communications, the signal-receiving apparatus 7 includes said 
cross-correlator that cross-correlates each received signal with said code 
sequence of the subscriber of which the transmitted signal is to be 
received. The output of the signal-receiving apparatus 7 forms the signal 
vector * and is connected to the signal processing apparatus 8. 

The signal processing apparatus 8 first generates the weight vector 
m by processing the signal vector x through the procedure presented in 
this invention. Then, the weight vector is used to generate the final 
output y of the antenna array system by performing the Euclidean inner 
product of the signal vector x and weight vector i.e., y = m H x, at 
each snapshot. The key part of the antenna array system shown in Fig. 1 
is the signal processing apparatus 8. It produces the weight vector w. at 
each snapshot that provides a nice beam pattern having its maximum gain 
along the direction of target signal source. 

Fig. 2 A is a flowchart that descibes the four steps of computing the 
weight vector m according to the first embodiment of the inventive 
method. The flowchart of Fig. 2A comprises the following steps. The 
initial estimation vector ^(0) and the initial gain vector w(0) are set 
with the initial signal vector ^(0) and the normalized estimation vector, 

i.e., j/<0) = ; *(0.) and _&?(0) = (| _^(o)n » and set the snapshot index n 
with 1, i.e., n-\ (step 21). The new signal vector &{n) is received, at 
the present snapshot (step 22). The estimation vector u(n) and the gain 
vector jm(ti) are updated, from the new signal vector by jy(w) = / 

M(n-l) + (l-f) An) x H (n) ju(n-i) and .jc(n) = . » respectively, 
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where the forgetting factor J is properly set with a value between 0 and 
1, i.e., 0£/<l (step 23). The final array output y(n) at the present 
snapshot is produced, by y{n) = w H (n) x(n) (step 24). Whether the 
signal processing is to be continued or not is checked (step 25). For 
continuing the procedure at the next snapshot, it goes back the step 22 to 
receive a new signal vector, with an incremented snapshot index «, i.e., 

« — nil (step 26). 

Fig. 2B is a flowchart that, descibes the six steps of computing the 
weight vector m according to the second embodiment of the inventive 
method. The main difference between the methods shown in Figs. 2A and 
2B is thai the target vector can be computed with repeated iterations at 
each snapshot in Fig. 2B, whereas the method shown in Fig. 2A 
generates the vector in a single iteration at each snapshot. The flowchart 
of Fig. 2B comprises the following steps. The initial estimation vector 

v(0) and the initial gain vector jcKO) are set with the initial signal 
vector ^(0) and the normalized initial estimation vector, respectively, i.e., 

y(0) =. *<0) and m>(0) = \\%(q)\\ > and the snapshot index n is. set with 

1, i.e., n = l (step 31). 

A new signal vector jf(«) is received, at the present snapshoht (step 
32). The estimation vector ji(n) and the gain vector w{n) are updated, 
(step 33), from the new signal vector z(n) by jitn) . = / y(n-\) +■ 
(1-y) xin) w(n-l) and win) = ? respectively, where the 

forgetting factor / is properly set with a value between 0 and 1, i.e., 0 <; 

/ <1. Whether or not more iteration is to be performed in the present 

snapshot is decided (step 33). Then, it goes to the following step 35 for 

more iteration in the present snapshot, or goes to the step 36 for no more 

iteration in the present snapshot. The estimation vector and weight vector 

are set (step 35), as ji(n-l) — ji(n) and m{n-\) «- " j^(«), 

respectively, for more iteration in the present snapshot. The final array 

output yin) at the present snapshot is generated (step 37), as y(n) = 

iv 1! (n) x{n) . And, whether the. signal processing is to be continued or 

net is check. For continuing the procedure at the next snapshot, it goes 

back the step 32 to receive a new signal vector, with an incremented 

snapshot index n, i.e., n <— w + 1 (step 38). 

Fig. 2C is a flowchart that descibes the four steps of computing the 

weight vector m according to the third embodiment of the inventive 
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method. The major difference of Fig. 2C compared to the previous two 
methods shown in Figs. 2A and Fig. 2B is that the autocorrelation matrix 
is computed in accordance with eq.(4) at each snapshot upon the reception 
of a new signal vector in the method introduced in Fig 2C. On the other 
hand, in the methods described in Figs. 2A and 2B the autocorrelation 
matix is not actually computed at each snapshot. Instead, the matrix is 
approximated with vector operations, as shown in eq.(13). Although the 
method shown in Fig. 2C requires more computation compared to that of 
Fig. 2 A or Fig. 2B, it can result in more accuracy in computing the 
weight vector and final array output. 

The flow chart of Fig. 2G comprises the following steps. The initial 
estimation vector j;(0) and the initial gain vector w(0) are set with the 
initial signal vector ^(0) and the normalized estimation vector, 

respectively, i.e., m(Q) = jc(0) and :w(0) = \\m(0)\\ ' anc * ^ e sna P s ^ ot 
index . n is set with 1, i.e., » = 1 .'"'(step 41). The autocorrelation matrix 
ll x {n) is updated, with a new signal vector x(n) received at the present 
snapshoht. by R x (n) = f R x (n— 1) + (1— /) x(n) _& H {n) where the initial 
matrix is set with the initial signal Vector by jg z (0) . = x(0) ^"(O) and 
the forgetting factor / is properly set with a value between 0 and 1, i.e., 
0:£/<l (step 42). The estimation vector ji(n) and gain vector win) are 
updated, from the new autocorrelation matrix by . u(n) - J2 x (n) win- 1) 

and m?( n) = || Jf( ^) j| > respectively (step 43). The final array output yin) at 

the present snapshot is generated, (step 44), as y(n) = jv"(n) jc(n) . And 
whether the signal processing is to be continued or not is checked (step 
450. For continuing the procedure at the next snapshot, it goes back the 
step 42 to receive a new signal vector, with an incremented snapshot 
index n, i.e.,. n n+l (step 46). 

Fig. 2D is a flowchart that descibes the six steps of computing the 
weight vector m according to the fourth embodiment . of the inventive 
method. The method described in Fig. 2D also computes the 
autocorrelation matrix at each snapshot, as shown in Fig. 2C. The major 
difference of the method shown in Fig. 2D compared to that of Fig. 2C is 
that repeated computation for the estimation vector and weight vector can 
be performed in the method described in Fig. 2D. Although the method 
shown in Fig, 2D requires more computation compared to the other 
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methods introduced in this invention, the method described in Fig. 2D can 
result in the best accuracy in computing the weight vector and final array 
output. 

The flowchart of Fig. 2D comprises the following steps. The initial 
estimation vector m(0) and initial gain vector _^(0) are set with the 
initial signal vector x(0) and the normalized initial estimation vector, 

respectively, i.e., v(0) = x(0) and _&K0) - ||Ji(0) II ' and the sna P shot 
index n is set with 1, i.e., w = l(step 51). The autocorrelation matrix 
_R x {n) is updated, with a new signal vector xin) received at the present 
snapshoht by R x (n) = f R x {n-\) + (1-/) x(n) *"<w), where the initial 
matrix is set with the initial signal vector by R x (0) = jz(0) z H (0) and 
the forgetting factor / is properly set with a value between 0 and 1, i.e., 
0:=c/<l (step 52). The estimation vector ji(n) and gain vector win) are 
updated, from the new autocorrelation matrix by. ju(n) = & x (n) Mt(n-l) 

and win) = ||^(^)|| » respectively (step 53). Whether or not more 

iteration is to be performed in the present snapshot is decided (step 54). 
Then, it goes to the following step 55 for more iteration in the present 
snapshot, or it goes to the step 56 for no more iteration in the present 
snapshot. The estimation vector and weight vector are set, as _£>(«— 1) — 
v(n) and -je(n-l) j&(n) , respectively, for more iteration in the 

present snapshot (step 55). 

The final array output y(n) at the present snapshot is generated, as 
y(n) = M H (n) ji(n) (step 56). And, whether the signal processing is to 
be continued or not is checked (step 57). For continuing the procedure at 
the next snapshot, it goes back the step 52 to receive a new signal 
vector, with an incremented snapshot index «, i.e., n *- w+l (step 58). 

Fig. 3 illustrates the performances of the antenna array system 
designed by the method presented in this invention. Fig. 3 A shows the 
number of clocks required to obtain the target weight vector by utilizing a 
microprocessor. In the simulations for the data of Fig. 3A, a digital signal 
processor, for example TMS320C30, has been selected arbitrarily. The 
required number of clocks in Fig. 3A is shown in comparison to another 
technique, namely a modified conjugate gradient method, that has been 
developed most recently. Fig. 3B shows the SIR and SNR performances 
of the antenna array system designed by the method presented in this 
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invention. 

As illustrated in the figure, the proposed system provides about 8 dB 
improvement in SIR. The SNR is, as mentioned earlier, improved linearly 
as the number of antenna elements in the array increases. From the 
performances shown in Fig 3, it is apparent that the antenna array 
system designed by the proposed method is excellent in terms of both 
performance and complexity. 

As mentioned earlier, the signal processing method provided in this 
invention gives the following advantages: first, the communication capacity 
is increased as much as the SIR are increased, and second, the 
communication quality is enhanced . as much as the SNR and the SIR is 
increased. 

The best feature of the proposed technique in this invention is that 
the required amount of computation to achieve all the merits is extremely 
small so that the proposed technique can be easily implemented with the 
norma] digital signal processor in real-time processing. The simplicity of 
the proposed method, as shown in Fig. 3, does not cause any sacrifice in 
accuracy compared to conventional methods. 

. Although the specific embodiments of the present invention have been 
disclosed and described, it is apparent that those who skilled in the art 
will appreciate that various modifications, additions and substitutions are 
possible, without departing from the scope and the spirit of the present 
invention as disclosed in the accompanying claims. Therefore, it should be 
understood that the present invention is not limited to the particular 
embodiment disclosed herein as the best mode contemplated for carrying 
out the present invention. 
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Claims 

1. A signal processing method for minimizing interference and for 
reducing effects of noise by controlling beam patterns of a 
telecommunication system having an antenna array, comprising the steps 
of 

(a) setting an initial estimation vector with an initial .. signal " vector,, 
setting an initial weight vector with normalized said estimation vector, and 
setting an initial snapshot index; 

(b) receiving a new signal vector at the present snapshoht; 

(c) updating said estimation vector and computing said weight vector 
by normalizing said estimation vector; and 

(d) producing a final array output at the present snapshot by 
computing an Euclidean inner product between the weight vector and the 
signal vector. 

2. The signal processing method according to claim 1, further 
comprising the step of (e) going back said step (b) to receive a new 
signal vector, with an incremented snapshot index, for continuing the 
procedure at a next snapshot. 

3. The signal processing method according to claim 1, wherein said 
step of (c) updating said weight vector comprises the substeps of: 

(c-1) multiplying each element of said estimation vector obtained at a 
last previous snapshot by a forgetting factor of which the value is 
predetermined to be in 0 to 1; 

(c-2) multiplying each element of complex conjugate of said signal 
vector received at the present snapshot by that of said weight vector 
computed during the last snapshot, and adding up all the results of the 
multiplications; 

(c-3) multiplying each element of said signal vector of the present 
snapshot by the result of said substep (c-2); 

(c-4) obtaining the updated estimation vector by adding each 
corresponding element of the results of said substeps (c-1) and .(c-3); and 

(c-5) determining said weight vector by dividing each element of said 
updated estimation vector, by a magnitude of said updated estimation 
vector. 
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4. The signal processing method according to claim 3, wherein said 
substep (c-1) comprises the following substeps of: 

(c-la) multiplying each element of said weight vector obtained at the 
last previous snapshot by said forgetting factor of which the value is 
predetermined to be iri 0 to 1;. and 

(c-lb) multiplying each element of the results of said substep (c-la) 
by the magnitude of said estimation vector obtained at the last previous 
snapshot. 

5. The signal processing method according to claim 1, wherein said 
weight vector is determined by multiplying a predetermined constant on 
each element of said estimation vector obtained in said step (c), in order 
to modify said weight vector without changing the beam -pattern 
characteristics of said estimation vector obtained in said step (c). 

6. The signal processing method according to claim 1, wherein said 
weight vector is determined by dividing each element of said estimation 
vector by the first element of said estimation vector obtained in said step 
(c), such that the signal received from a reference antenna element is 
always multiplied by a real -valued gain. 

7. The signal processing method according to claim 6, wherein said 
weight vector is determined by normalizing the result of the division, such 
that the magnitude of the final resultant weight vector is 1. 

8. The signal processing method according to claim 7, wherein said 
reference antenna element is determined by an antenna element of which 
the phase of said signal is the latest of all said antenna elements in said 
antenna array at the present snapshot. 

9. The signal processing apparatus according to claim 7, wherein said 
reference antenna element is determined by said antenna element of which 
the physical distance from a signal source to be communicated with at 
the present snapshot is farthest compared to the other antenna elements 
in said antenna array. 
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10. The signal processing method according to claim 1, further 
comprising the step of (e) going back to said step (c), with an updated 
weight vector substituting for the previous weight vector, for more 
iteration in the present snapshot. 

11. The signal processing method according to claim 2, further 
comprising the step of (f) going back to said step (c), with an updated 
weight vector substituting for the previous weight vector, for more 
iteration in the present snapshot. 

12. A signal processing method for minimizing interference and for 
reducing effects of noise by controlling beam patterns of a 
telecommunication system having an antenna array, comprising the steps 
of 

(a) setting an initial estimation vector with a signal vector received at 
initial snapshot, setting an initial weight vector with a normalized 
estimation vector, and setting an initial snapshot index; 

(b) updating an autocorrelation matrix with a new signal vector 
received at a present snapshot, where an initial matrix, is set with the 
initial signal vector and a forgetting factor is set with a value between 0 
and 1 '» 

(c) updating the estimation vector and the weight vector from the 
new autocorrelation matrix, respectively', and 

(d) producing a final array output at the present snapshot by 
computing an Euclidean inner product between the weight vector and the 
signal vector. 

13. The signal processing method according to claim 12, further 
comprising the step of (e) going back said step (b) to receive a new 
signal vector, with an incremented snapshot index, for continuing the 
procedure at a next snapshot. 

14. The signal processing method according to claim 12, wherein said 
weight vector is determined by multiplying a predetermined constant on 
each element of said estimation vector obtained in said step (c), in order 
to modify said weight vector without changing the beam -pattern 
characteristics of said estimation vector obtained in said step (c). 
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15. The signal processing method according to claim 12, wherein said 
weight vector is determined by dividing each element of said estimation 
vector by the first element of said estimation vector obtained in said step 
(c), . such that the signal received from a reference antenna element is 
always multiplied by a real-valued gain. 

16. The signal processing method according to claim 15, wherein said 
weight vector is determined by normalizing the result of the division, such 
that the magnitude of the final resultant weight vector is 1. 

17. The signal processing method according to claim 16, wherein said 
reference antenna element is determined by an antenna element of which 
th? phase of said signal is the latest of all said antenna elements in said 
antenna array at the present snapshot. 

18. The signal processing apparatus according to claim 16, wherein 
said reference antenna element is determined by said antenna element of 
which the physical distance from a signal source to be communicated with 
at the present snapshot is farthest compared to the other antenna 
elements in said antenna array. 

19. The signal processing method according to claim 12, further 
comprising the step of (f) going back to said step (c), with an updated 
weight vector substituting for the previous weight vector, for more 
iteration in the present snapshot. 

20. The signal processing method according to claim 13, further 
comprising the step of '(g) going back to said step .(c), with an updated 
weight vector substituting for the previous weight vector, for more 
iteration in the present snapshot. 
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